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Abstract 

We calculate the probability of the decay of external inhomogeneous electro- 
magnetic fields to neutral pseudoscalar particles that have a coupling to two 
photons. We also point out that our estimate for axion emission in a previous 
paper was incorrect. 
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I. INTRODUCTION 



The Schwinger mechanism is a non-perturbative process by which an infinite number of 
zero frequency photons can decay into electron-positron pairs 0. In this paper we show 
that this mechanism can be generalized to study the production of other kinds of light 
particles from intense electromagnetic (EM) fields. The light particle that we consider is a 
pseudoscalar (PS) having a coupling to two photons. 

In Section II we derive the formula for the decay of classical background fields into PS 
particles. This is achieved by integrating out the particle fields from the total Lagrangian 
to obtain the effective action of the classical background fields. The imaginary part of the 
effective Lagrangian is related to the probability of decay of classical background fields into 
particles. In Section III, we derive, from the usual coupling of the PS to two photons, the 
specific interaction Lagrangian that should be used in the general formalism of Section II in 
order to account for vacuum decay into PS. For static EM fields, we show that a necessary 
condition is that the fields are inhomogeneous. In Sections IV, V, and VI, we explicitly 
calculate the PS production in a variety of situations. Specifically we consider a dipole 
magnetic field, a cylindrical capacitor, and a spherical capacitor. A final section is devoted 
to the conclusions. 



II. DECAY OF CLASSICAL BACKGROUND FIELDS INTO PARTICLES 



We start with the action for the pseudoscalar (mass m) coupled to the background E 
and B fields of the general form 



S[<P,E,B] 



d A x -<f)(x) —d 2 — m 2 + f(x) <j)(x) 



(1) 



where f(x) is some scalar function of E and B fields. From ([!]) we obtain the effective action 
for the background E and B fields formally as 



V<f> e l 



■,E,B1 



(2) 



The effective Lagrangian for the E and B fields can be related to the Green's function of 
in external E and B fields as follows. Differentiate @) by m 2 

.dS eff [B,B) f V<j)<p 2 e lS ^ E ^ 



dm 2 



fV(f) e i5 ^' E < B ] 



-I I d 4 iG(ix;E,B) 



dp 



G(p;E,B) 



(3) 



The effective Lagrangian of the background fields is therefore formally given by the expres- 
sion 

d 4 p 



£ e //[E,B] 



dm 



(2tt)< 



G(p;E,B) 



(4) 
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The probability of external E and B fields to decay into quanta of 
imaginary part of £ e ff as follows 



is related to the 



P = l-(0 
= 1 — exp 



»S e //[E,B] 







2 Im / d 3 x dt C eff [ 



E,B1 



(5) 



In the case that this probability is small, we can write the probability density w (per unit 
volume and unit time) approximately as 



w = 2 Im£ e// [E,B] 



(6) 



We now give the general procedure for obtaining the effective action of the background 
fields by calculating the Green's function of in background E and B fields following the 
method of Duff and Brown ||. 

The effective Lagrangian can be calculated by this method if the background fields con- 
tained in f(x) in the interaction Lagrangian 



Cj(x) = l -f{x)<p\x) 



(7) 



can be expanded in a Taylor series near some reference point x. Expanding f(x) near x = x, 



f(x) = a(x) + (3^{x){x - xY + % v (x)(x - xY(x - x) v + ... 



(8) 



a (x) = f(x), Pn(x) 



'or 

dx^ 



2 If d df 



2 \ dx^ dx h 



The equation for the Green's function for the <fi field is given by 



dl + m 2 - a - P^x - xY - 7^,(2 - x)^(x - x) u G(x, x) = 5 4 (x - x) 



In momentum space 



fx — xY — * —i 



_d_ 

dp. 



and the equation for the Green's function in momentum space is 



22 -n 9 2 dd 

-p +m -a + iPn- h7 - 

dp^ " op^dpv 

We choose as an ansatz for the solution G(p) the form 



G(p) = 1 



G(p) =1 ds e -is(m 2 -ie)^ Pf ,A^p,+B^+C 
JO 



(9) 



(10) 



(11) 



(12) 



where A(s), B(s) and C(s) are to be determined. They must satisfy the boundary condition 
in the case of vanishing external fields, i.e. when a, (3, 7 —>■ 
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M lv -> s g"", 
and in this limit we should obtain 



B» -> 0, 



/•DO .9.9 

G(p) = i / ds e - lsm +lsp 
Jo 



m 2 — p 2 



(13) 



(14) 



i.e., the free particle Green's function. 

To solve for A, B, and C we insert ansatz (O) in (TTTT) . We have 



/■oo 

i / <is 
Jo 



- p 2 + m 2 -a + i/3- (2iA ■ p + B) + (2ip ■ A + B) ■ 7 2 ■ (2zA ■ p + B) 



+ 2z tr(7 2 • A) exp {-ism 2 + zp-A-p + B- p + c} = l 
Equation (|i~5"D has the general form 

/ ds g(s) e~ h{s) = 1 

whose solution is 

dh(s) 



ds 



(15) 



(16) 



(17) 



with h(0) = and h(oo) = 00. Using the form of the solution ( [l~7f) for equation (|15| 
- p 2 + m 2 - a + i(3 ■ (2iA ■ p + B) + (2ip ■ A + B) ■ -y 2 ■ (2iA ■ p + B) 
+ 2ztr( 7 2 ■ A) 



<9A 



9B 



«m - «p • — • p - — • p 

OS OS 



dC 
ds 



(18) 



and comparing equal powers of p on both sides we get the following linear differential equa- 
tions for A, B, and C, 

<9A , , , 

— = 1 + 4A • 7 2 • A 

as 



OB 

— = 2iA ■ (3 + 4A • 7 2 ■ B 

OS 



dC 
ds 



ia + (3 ■ B - iB ■ ~/ 2 ■ B + 2 tr(7 2 • A) 



(19) 



The solutions of these equation which satisfy the boundary conditions (13|) are given by 



A = - 7 1 • tan(27s) 



B 



7" 2 • [1 -sec(27s)] • (3 



1 i 

C = ias - - tr [In cos(27s)] + - (3 ■ 7~ 3 • [tan(27s) - 27s] • (3 
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(20) 
(21) 
(22) 



These A, B and C determine G(p) when substituted in ([121) . The effective Lagrangian is 
obtained by substituting this G(p) in (f|) and carrying out the integration over m 2 , 



i f 00 ds f cTp r . 2 _ 



d 4 p 



(2vr 



The Gaussian integral may be evaluated using 



d 4 p exp {ip • A ■ p + B • p} = — in 2 (detA) 2 exp 



- B • A" 1 • B 
4 



(23) 



(24) 



where detA is the determinant of the matrix A%. Using now (P0| - P2| ), we have 



ef! 



327T 2 Jo S 



ds 



-is(m 2 —a) 



det 



'Jrysj 
sin 27s 



(25) 



where 



l( S ) 
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f3 ■ 7 3 ■ [tan(7s) - 7s] • (3 



(26) 



The coefficients of the Taylor expansion of the background fields (§) determine the effective 
action on integrating out the quantum field (p. In particular, an imaginary part of £ e // m &y 
be non-zero depending on the signs of the eigenvalues of the r ) 2 matrix. When this occurrs, 
we have a non-zero probability (^) that the external EM fields decay in PS particles. 

To the effective Lagrangian in ( p5[) we should add subtractions to render it finite at 
s = 0. When this is done, we have that in the limit (3 —>■ 0, j 2 — > 0, the effective Lagrangian 
£ e ff — ► 0, as it should be. In the Appendix A we illustrate the method for the familiar case 
of production of charged scalar fields in a constant electric field. 

The formulae ( ptj| - p2| p5|) differ by some signs and factors of i from the solutions displayed 
in reference 0. There, we presented the formulae for the case that f(x) had only spatial 
variation and therefore (3 and 7 2 had only i = 1,2,3 indices. In [H] we used the metric 
(+, +, +) while in the present paper we consider spatial as well temporal variation and 
use the metric (+, — , — , — ). This introduces some changes in intermediate formulae but of 
course the final results we get in the present paper are identical with the final results we got 
in 0. 



III. EFFECTIVE EM-FIELDS - PS PAIR INTERACTIONS 

In this Section, we show that a coupling of a pseudoscalar to two photons induces an 
interaction Li that may lead to PS production in a background of EM fields. 
The generic pseudoscalar- two-photon interaction (see fig.l) can be written as 

= l -g<Pe^F^F pu (27) 

We should mention that in the special case where the PS is an axion the coupling g is related 
to the mass m of the axion. 

We need to evaluate the loop diagram of the type shown in fig. 2 with infinite number of 
zero-frequency photon external legs. The imaginary part of this diagram gives the probability 
for the decay of the external electromagnetic field. 
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To calculate this diagram, we first evaluate the process <pA — > (f)A, where A is an external 
photon. We use from fl27D in momentum space, 

- A g^°k,A u F pa (28) 
The two-photon two-PS interaction is then obtained contracting the internal photon legs, 

4 (j#)^ M ^^^^ Vff '(-V)^' (29) 

The factor of 4 in equation (|29| ) is for the four possible ways of joining the photon legs. Due 
to the presence of the k 2 term in the denominator, the effective coupling (|29"D is non-local. 
However, when we calculate the effective action for the external EM field the momentum k 
is integrated over. One can therefore make use of the identity 

J d'k kfy g{k 2 ) = J d 4 k 9 -^- g{k 2 ) (30) 

to simplify ( p9j) . Thus, we can reduce the effective two PS-two photon interaction to a local 
interaction vertex. Back in configuration space, it is given by 

Cj = -Ig^F^ = ^ 2 2 (E 2 - B 2 ) (31) 

(see fig. 3). 

With the interaction Lagrangian ( |31~D we can go back to the formalism of Section II and 
calculate the probability density. We can readily identify f(x) in ([/]), 

f(x) = g 2 (B 2 -B 2 ) (32) 

In order to have a non trivial C e ff, one needs non-zero second derivatives of the EM fields 



as they appear in expression ([32]) . As we said in Section I, depending on the sign of the 
corresponding r y 2 matrix we may have PS production. We illustrate it in some simple physical 
situations in the following sections. 

IV. PRODUCTION OF PSEUDOSCALARS IN DIPOLE MAGNETIC FIELDS 

In a static dipole magnetic field the PS-pair -EM interaction is given by 

=4 92 ( B -» !£ ¥^ ! )* 2 (33) 

where B is the field strength at a point r = (0, 0, z ) on the z-axis. We have now 

m = -J "{bI^-^A (34) 



6 



Expanding B 2 (r) near the point fo 



C 



-d 2 -m 2 + /(f ) + 



df_ 

dxi 



+ 



1 d 2 f 



2 dxidxj 



{Xi - Xio) 

f=f 

(Xi - X i0 )(Xj - X jQ ) (f) + 



r=ro 



we find that the coefficients of the Taylor expansion are given by 

a = f(ro) = ~g 2 Bl = a m 



(0,0, e^oV 



r=r 



and 



(7 2 



1 d 2 f 



i i.i 



2 dxidxj 



r=r 



2 


(5 







( aL 
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:). 





a 2 m 




v0 





-28 / 


V o 


-b 2 m 



(35) 
(36) 

(37) 
(38) 



(39) 



Therefore we find using the notation and formalism of Section II that the effective action 
on integrating out the PS field is given by 



C 



eff 



32vr 2 Jo s 



00 ds 



-is(rn 2 -a m ) ^ a mS 



2b m s 



sinh 2a m s y sin 2b m s 



(40) 



with 



lm{s) = Xm{b m s - tan b m s) 



(41) 



9g 4 B* 1 3 „ 



The imaginary part of the expression ( f4"0"D can be performed by enclosing the simple poles 
at s = —inir(2a m )~ l , n = 0, 1, with a contour from below. We get 

1 3 i OO 

Im C eff = —alblY. (-l^C^e-™^™ (43) 
8tt2 n=1 



(j(m) 



n 2 



sinn n — it 



(44) 



with 



15 / m 



m 2 - a m 2 \ gB 



z gB 



(45) 
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and 



b m 7l b m TT 

n tannn 

dm. 2 &rn. 2 



(46) 



In (|44D and (f46|) we can put b m /a m = y28/5. The main contribution to the above integral 
comes from the n = 1 term and we find that w is given by the expression 

w = .036^fV- m - e~ n/2 ^ (47) 

The probability w of field decay is extremely suppressed for realistic parameters. To 
illustrate this, let us choose a mass m and a coupling g consistent with the axion window: 

m ~ 10" 3 eV 

g ~ 1(T 13 GeV" 1 (48) 



Also, let us choose 



Bq = 1 Tesla 

z = 10 cm (49) 



1(T 19 (50) 



With these values, we get 



and since rj m appears in the exponential in (|47j), we see that the probability w is indeed 
extremely suppressed. For a dipole magnetic field to be unstable and decay into axions, one 
needs rj m ^ 1, but this would correspond either to unrealistic values for the external field 



parameters (f49D or to excluded values for the axion mass and coupling (f48D . For non-axion 
models, g and m are not related (still there are restrictions on these parameters, see ref. M). 
One could have r] m ~ 1 by tuning g and m. Imposing that the field fl49|) does not decay into 
pseudoscalars leads to the constraint 



(51) 



m \ > 9 



10- 12 eV/ 10- 13 GeV- 



V. PRODUCTION OF PSEUDOSCALARS IN A CYLINDRICAL CAPACITOR 

The modulus of the electric field inside a cylindrical capacitor whose axis lies along the 
z-axis depends only on p = (x 2 + y 2 ) 2 , 

EM = I'- (52) 



with A the linear electric charge density. 
The bilinear interaction term (|3l|) is 



S 



C, = ij 2 B 2 (p)^(i) 

= 5* (7) *"<*> 
where g c = Xg/27T. The corresponding function /(x) is 

f(p) = 9l (^) 

Expanding the fields near some reference point (xq, Uo,Zo) with p = O^o + 2/o) 5 



(53) 



(54) 



-d 2 -m 2 + f(p ) + 



df 
dp 



(p - Po) + 



Id 2 / 



P=PO 



2 dp 2 



(P-Pof 



P=P0 



+ .. 



(55) 



It can be written as in (H) with 



« = /(Po) = «c 



(56) 



and 



(7 2 ) 



f 

(p)i = —(x ,y ) 
Po 



2pl \-x y x 2 J 2p 2 Q \x y y 2 J 



(57) 



(58) 



where primes denote derivatives with respect to p taken at p = p . In the above formulae, 
the spatial indices run over 1,2. 

Introducing the explicit form of /, we get 



Or 



pI 



(59) 



2gL , 
■—{xo,yo) 

Po 



(60) 



and the r y 2 matrix (|58|) reads 



.2\ 9cf-y 2 + ^x 2 Ax y 
1 )ij p% \ 4x y -A + H 



(61) 



Next, we diagonalise 7 2 by rotating the coordinates with an orthogonal matrix. For 
example, we can use 



1 '' %o Vo 



(Oh = - , 

Po v -yo x o 



(62) 
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In diagonal form we have 



We need /3 in the diagonal basis given by ftp = O • (3. We get 

{P D )i=(-2&Po 3 ,0) (64) 
The expression for C e ff ( p5|) finally reads 



eff 32Wo s 3 Vsinh2a cS Vsin26 cS 1 j 



where 



l c (s) = A c (a c s — tanha c s) (66) 
9c 
3V3 



Ac = gtPoV = (67) 



We are not able to perform the integration in (|65| ) by the procedure of extending s to 
the complex plane, as we have done in section IV. The reason is the presence of essential 
singularities contained in l c (s). (For a discussion of the implications of essential singularities 
in the context of QED pair production calculations at finite temperature, see ref.( |J.) 

We shall calculate the integral (|6~5|) numerically. We make the change of variables 



so that 



x = (m — a c — X c a c ) s (68) 



w = 21m C ef f = 1 — -zl c (69) 

IQir 2 (m 2 — a c — \ c a c ) 2 



where 



00 



g?x . t ± N / / 2rf c x 
x 3 



„3 Sln ^c; 



2r/ c a;/V / 3 , 2 



+ -y - 1 ( 7 °) 



sinh277 c x \ sin(r7 c x/\/3) 9 
We have introduced the necessary subtractions and defined 

Vc = 2 ^ r (71) 
771 — a c — A c a c 

and 

0c = x + A c tanh ?7 c x (72) 

As it stands, 7 C depends on the two adimensional parameters A c and r] c , which reflect its 
dependence on both the strength of the interaction and on the mass of the scalar particles. 
In order to explore numerically this two-dimensional space we should focus on those regions 
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for which the results make physical sense. This is even more so because a blind computation 
of the integral leads easily to wild fluctuations due to the violent oscillations of the integrand 
for large domains of parameter space. 

Since the electric field should be inhomogeneous on scales of the order of a particle's 
Compton wavelength, we have 




> 1 



In our case, this means 



p m < 1 



(73) 



(74) 



Moreover, we should restrict our survey to subcritical conditions, i.e. conditions such that 
we find ourselves below the point where catastrophical pair-production starts in the earnest 
and vacuum breakdown occurs. The system will be subcritical whenever the field energy 
stored in a volume of size m -3 , and which is converted into a PS pair, is at most of the order 
of twice the scalar particle rest mass. A crude back of an envelope estimate gives, 



Xl/mV < 0(1) 



(75) 



Using both restrictions as a rough guide, we perform the numerical integration of J c as 
a function of r]~ l for i] c and A c in the ballpark of the values required by (|74|) and (75). The 
results for two values of A c are displayed in figs. 4 and 5. The curves are accurately fit by an 
analytic expression of the form, 



Arjle 



-k/rj c 



(76) 



with A and k depending on A c . (A and k are positive.) Equation (^) has a form that 
closely resembles the classical Schwinger result and characterizes a typical non-perturbative 
process. 

The electric field would break down into pseudoscalars if the exponent in ( [75] ) becomes 
small. Again we choose the values consistent with the axion window (^) and for the field 
parameters we take 



E ~ 10 4 V/m 
p ~ 0.1 cm 



(77) 



This leads to 



Vc ~ 1(T 22 ( 78) 

which implies a large negative value of the exponent in (|76|) that suppresses field decay. 



VI. PRODUCTION OF PSEUDOSCALARS IN A SPHERICAL CAPACITOR 

The modulus of the electric field inside a spherical capacitor depends only on r = \ f\ 
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where Q is the electric charge. 

The bilinear interaction term is then 



Cj = -g 2 E 2 (r)<p 2 (x) 



= ~ 2 9t ) n*) 

where g s = Qg/An. 

The corresponding function f(x) depends only on r, 



M = 9 2 s Qr 4 



Expanding the fields near some reference point with (xq, 2/o, ^o) with modulus ro 



C 



-d 2 - m 2 + f(r ) + 



Or 



r=r 



[r - r ) 



r=r 



+ 



In Cartesian coordinates it can be written as in (|8|) with 

a = f(ro) = a s 



f 

(p)i = — (x ,y ,z ) 
r 



and 



(7 2 



1 1.1 



JL 

2rl 



—XqZq 



x + z 



— XqZq 

-yoz 
A + vo 



+ 



r_ 

2rl 



I x, 



XqVq 

\XqZq 



where primes denote derivatives with respect to r taken at r 
f(x), we get 



XoVo x z 

VoZo 4 
r . Introducing the 



a-. 



and 



and 



h 2 ) 



2^ 

.',-.8 



6x y 
\ &x z 



Ag 2 

—£-(x ,yo,z ) 



—x + 5y — z 



6x z 
fyoZo 
~ x o ~ Vo + ^z 
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The 7 2 matrix ([88]) can be diagonalised by rotating the coordinates with an orthogonal 
matrix. For example, we can use 

I ( -ZqTq x r 
(O)y = —r -x y d 2 Q -y z I (89) 
r ° 1 x d y d z d 



iiDa = % I 1 I = 














-b 2 s 











a* 



(90) 



where do = \Jx1 + z\. 

In diagonal form we have 

-1 0' 
) -1 
5 

One must also use (3 in the diagonal basis given by [3d = O ■ (3. We get 

(J3 D )i= (o,0,-4<7 c V) (91) 
The expression for C e ff (]25|) is given in this case by 



c = _J_rds e ^ as) 2a s s _2M e «, w (Q2) 
/; 327r 2 7o s 3 Vsinh2a s s sin26,s v ; 



where 



Z s (s) = A s (a s s — tanha s s) (93) 
4ff s 4 1 2 & 
5 V^O r 

As in the precedent section, we cannot perform the integration in ( |92"D by extending s 
to the complex plane since again there are essential singularities and we do a numerical 
integration. The procedure is very similar. It is convenient to change variables, 

x = (m 2 — a s — X s a s ) s (95) 



so that 



where 



°°dx . f I 2 Vs x 2 Vs x/V5 1 2 \ 



sinh2r7 s x sin^x/v^) 5 

We have introduced the necessary substractions and defined 

a s 

V s = ~ \ 



(98) 
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and 



4> s = x + \ s tanh r] s x (99) 

Here we follow a similar strategy as before to pinpoint the relevant parameter space. We 
get similar restrictions: 

r m < 1 and \ 2 Jr 2 m 2 < 0(1) (100) 

In figs. 6 and 7 we display I s as a function of rjj 1 , for a couple of values of A s . Again, we see 
that the approximate behaviour is that of a decreasing exponential, 

I s = Arf s e~ kl ^ (101) 

with positive constants A and k (that depend on A s ). 

When one considers realistic values for the field and PS parameters it turns out that the 
probability of field breakdown is extremely suppressed as it happens in the cases that have 
been analyzed in the precedent sections, namely the dipole magnetic field and the cylindrical 
capacitor. 



VII. CONCLUSIONS AND FINAL REMARKS 

In the presence of strong external fields, the physical vacuum breaks down because 
particle-antiparticle pairs are being pumped out of it at the expense of field energy. The case 
of a strong uniform electric field spontaneously creating electron-positron pairs is the best 
known (QED) example for this phenomenon. Such process is of a non-perturbative nature 
and the QED case has been solved exactly by Schwinger and others [[TJ. Their solution 
however does not include the backreaction on the external field exerted by the presence of 
the produced e + e~ pairs. Clearly, creation of pairs requires the supply of mass energy and 
kinetic energy which must be furnished by the external field. A balanced energy budget is 
therefore only possible through a corresponding reduction of the energy stored in the field. 
Because electrons and positrons carry charge they will fly to the external sources of the field 
and thus the field (and hence its energy) will diminish. So, unless from the outside the field 
is restored, the pair production process cannot be indefinitely sustained. If nothing is done 
from the outside a catastrophic breakdown of the initially strong (critical) electric field will 
inevitably follow. 

In the present paper we dealt with pseudoscalar particles. Pseudoscalars are fundamental 
ingredients of many completions of Particle Physics models. Examples run from axions to 
superlight partners of gravitinos. In the previous sections we have derived the probability for 
pair production of PS in electric and magnetic fields. Contrary to the QED case mentioned 
above, constant fields do not cause the disruption of the vacuum. Field gradients are neces- 
sary for the phenomenon to occur. Hence, we studied PS pair production in inhomogeneous 
fields. We have calculated the probability in a general case and based our computation on 
an effective action formalism formulated by Brown and Duff. We then have applied the gen- 
eral formulae to a few specific cases: PS production in a magnetic dipole field and between 
the plates of a charged capacitor (either cylindrical or spherical). Again, backreaction was 
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ignored and therefore adequate boundary conditions were implicitly assumed that take into 
account the fact that pairwise creation of PS requires field energy to be depleted. 

In the three cases studied, we found that our probability shows the non-perturbative 
behaviour w ~ exp(— const x m? / g) expected for subcritical fields. Finally, we should point 
out that in a previous paper || we erroneously estimated axion emission in the Coulomb 
field of an atomic nucleus. This result is incorrect because we overlooked the question of 
appropriate boundary conditions that guarantee energy conservation and which are clearly 
not met in this microscopic system. 
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APPENDIX A: DECAY OF A CONSTANT ELECTRIC FIELD INTO CHARGED 

SCALARS 



We start with the equation for the Green's function G(p) 

-(p - eA) 2 + m 2 ] G(p) = 1 
-p 2 + m 2 + e(A% + ppA") - e 2 A 2 ] G{p) = 1 

We assume constant E and B fields. The vector potential can be choosen as 

A =--F x^-F 9 



2 ^dp„ 



(Al) 



(A2) 



When inserted in ( |Al|) , one gets an equation of the form given in (|TT|), except for a term 



F^^-G(p) 

Op v 

that leads to an expression containing 

F^ PflPa A^ 



(A3) 



(A4) 



The antisymmetry of F and the fact that F and A commute makes this term vanish. Then 
our equation is 



e- d 2 



- p2 + m2 + ^ F " d p,Op p 



G{p) = 1 



(A5) 



and with our definitions in (181), (3 = and 
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lup ^ 1 liv 1 p 

e 2 

= --F 2 

4 v p 

Let us work out the special case of a constant E field. We have 



IE 2 



(7 2 ) 



vp 



e 

T 



o 



V o 









\ 




-E 2 ) 



(A6) 



(A7) 



(We have chosen the ^-direction as the direction of E). The eigenvalues of j are negative 
so 



det 



2js x 
sin 27s 



eEs 



sinh eEs 



(A8) 



and 



C 



1 



eff 



327T 2 JO s 



00 ds 



eEs 



sinh eEs 



(A9) 



To this expression for £ e // one should add a subtraction to make it finite at s = 0. When 
this is done, £ e // — > when eE — > 0. 

The probability of scalar production can now be calculated using (^|). The integral can 
be calculated by contour integration by closing the real axis with a contour on the negative 
imaginary plane. This contour encloses poles at s = —inn/eE which contribute to the 
integral. The final result for the constant electric field decay probability density is 



w 



n/:,' 2 g 



-1 



m+l 



2tt 2 



n=l 



n- 



exp 



eE 



(A10) 



which coincides with the well-known formula found in textbooks 
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FIG. 4. I c as a function of l/rj c for the value A c = 0.1. Dotted line: I c obtained by numerical 



integration of (70). Full line: I c as given by (ffq) with A = 0.28 and k = 1.026 
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FIG. 5. Same as fig.4 for A c = 0.4 and A = 0.43 and k = 2.1 
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FIG. 6. I s as a function of \/r} s for the value A s = 0.2. Dotted line: I s obtained by numerical 
integration of (p7|). Full line: I s given by (101) with A = 0.30 and k = 1.46 
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FIG. 7. Same as fig.6 for A s = 0.5 and A = 0.35 and k = 1.92 
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